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We present a geometric and dynamical approach to the micro-canonical ensemble of classical 
Hamiltonian systems. We generalize the arguments in and show that the energy-derivative of 
a micro-canonical average is itself micro-canonically observable. In particular, temperature, spe- 
cific heat and higher order derivatives of the entropy can be observed dynamically. We give per- 
turbative, asymptotic formulas by which the canonical ensemble itself can be reconstructed from 
micro-canonical measurements only. In a purely micro-canonical approach we rederive formulas by 
Lebowitz et al relating e.g. specific heat to fluctuations in the kinetic energy. We show that un- 
der natural assumptions on the fluctuations in the kinetic energy the micro-canonical temperature 
is asymptotically equivalent to the standard canonical definition using the kinetic energy. 
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I. INTRODUCTION 



For an isolated classical Hamiltonian system the ergodic hypothesis asserts that the time average of an observable 
CNJ ■ along almost any trajectory may be replaced by a space average over a suitable subset of the phase space, typically 
the energy surface. Such an average is denoted the micro-canonical ensemble average or for short, /z-average. The 
thermodynamic variables in this ensemble are the first integrals as well as extensive quantities like volume and particle 
numbers (cf e.g. Abraham and Marsden Chapter 3.7]). In the so-called thermodynamic limit of many systems 
, coupled weakly one expects in equilibrium each individual system to behave according to the canonical or Gibbs 
ensemble. In this ensemble the free parameters J/j are the variables conjugated to (some of) the first integrals. Quite 
f~^. | simple statistical principles argue in favor of such an approach but rigorous results are sparse || . Whereas the theory 
0^ • of the canonical ensemble has been elaborated to almost perfection, making it a corner stone in modern physics, our 
present understanding of the /^-ensemble and the equivalence of the two ensembles is remarkably incomplete. On the 
other hand modern computer technology makes it ever more important to understand and give precise interpretations 
' for dynamical measurements in the /x-ensemble j5j . 

From a geometrical point of view the /i-ensemble is an average of smooth quantities over a (genericly smooth) 
sub-manifold of phase space, fixed by the first integrals. For the sake of clarity we restrict ourselves to the case 
where the energy is the only extensive parameter. But even in this simple picture it is not a priori clear how to 
O ■ describe the important role played by the variable conjugate to the energy, i.e. the inverse temperature. When the 
Hamiltonian contains one or more separate terms of the form momentum squared the canonical average of such a term 
• »-H , yields precisely the (canonical) temperature. The time-average of such a term is therefore often used as a measure of 
the physical temperature of the system (cf. Example 3.7.27] or ||). This approach, however, presumes both the 
ergodic hypothesis and the equivalence of ensembles. In [ fiof we used simple geometrical arguments to show that the 
//-temperature is in fact measurable in the /i-ensemble itself. In particular, assuming ergodicity only, we constructed 
explicitly an observable whose average yields the /i-temperature. In the present article we shall show that these 
arguments in fact carries over to a much wider range of /Lt-observables. 

Below we establish (Theorem 1) two fundamental identities that will allow us to measure any energy-derivative of 
a /^-average within the /x-ensemble itself. Thus not only the temperature but also the specific heat and any higher 
order derivatives of the /i-canonical entropy can be observed dynamically. 

We give two main examples to illustrate these principles. In the first we use a perturbative expansion to show that 
by measuring all energy derivatives in the micro-canonical ensemble we can in principle reconstruct the canonical 
ensemble, cf. equations (|2|) and (|2~i|). This approach is, however, based upon a Gaussian expansion which itself relies 
on analyticity. On general grounds one would expect such results to be at best asymptotic, e.g. in the number of 
degrees of freedom. 

In the second example we consider the thermodynamics of particles in a box. Lebowitz et al || established relations 
between e.g. the specific heat and fluctuations in the kinetic energy in the micro-canonical ensemble. We rederive 
these relations in a purely micro-canonical approach. We also show that when both the kinetic energy and the square 
of fluctuations in the kinetic energy are of order N our micro-canonical derivation of temperature is asymptotically 
equivalent to the standard canonical definition of temperature which uses the average kinetic energy. It is of interest 
to study deviations in the two approaches for systems out of equilibrium. 
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II. ENERGY DERIVATIVES 



The phase space is an n — 2A-dimensional symplectic manifold (A4,uj) where u is a non-degenerate symplectic 
two-form and the Liouville volume form m = /\ n uj is nowhere vanishing. The standard example is (a subset of) 
Euclidean space M. = R 2N and m = Lebesgue measure (cf. Example B below). A Hamiltonian function H : M. — > R 
generates jl],|| then a vector field IdH and a flow g\ dH which preserves the Liouville measure as well as the energy 
surface, He = {£ G A4 : = E}. It follows ( fjL chapter 3.7]) that the flow also preserves the restricted Liouville 

measure, formally given by : 

H E = m5{H -E). (1) 

The /^-measure of an observable cf> is given by he(<I>) — J m 6(H — E) <f> and is a function of the energy. The 
/i-ensemble is the corresponding probability distribution, yielding the /^-average of </> over the energy surface : 

(0;£)= ME (0)/ MJ5 (1). (2) 

If the Hamiltonian flow is ergodic with respect to the restricted Liouville measure then by Birkhoff's theorem this 
/Lt-average equals the time average of <f> for almost any initial point on the energy surface Tie- We say that ((/>; E) is 
measurable in the /i-ensemble. 

Of particular importance is the /x-entropy and the associated (inverse) temperature : 

S^-Wh m 1 dS(E) 

s(£)-io g Mi), ( 3 ) 

The energy derivative and the /i-averages can now be related through the following : 
Theorem : 

Assume that He is a regular energy surface of the Hamiltonian function H and that X is a vector field defined 
in a neighborhood of He satisfying : 

dH(X) = 1. (4) 
Then the ^-measure and the ^-average of a observable <fi verify the identities : 

^e(0) = Mdiv (0A)), (5) 

±{^E) = (dxv(4>X);E)-^. (6) 



In coordinates, X • V = Xj^7, m = p dx\ A ... A dx ni one has the explicit formulae : 

1/9 f) FT 

div(^X) = -^— (rfXi), dH(X) = Y^Xi— = l. (7) 

It is sufficient to verify the identity (|^) since the other follows from (H) and 

dE K9 ^> M£ (l) Ml) ' 1 j 

According to Khinchin the /i-canonical measure is proportional to dE/||VJ?|| where d-H is the area element on 
the energy surface. Thus for an observable <f) we have 







/ m<f) — 


[ du [ 


IH<E 


l-oo JH 
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In the language of differential geometry the /i-canonical measure can be expressed as a differential (2n — l)-form : 

fi = i x m, (10) 

where A is any vector field for which dH{X) = 1. Provided H has no critical points on the energy surface one can 
always find such a vector field in a neighborhood of that energy surface. Although /i is not unique, its restriction to an 
energy surface is unique and equivalent to the weighted area element in Khinchin's formula. The (exterior) derivative 
of n<p {<f> being an observable) is given by : 

d{jjb(f>) — d(ixm4>) = dfa^xrn) = rn div((f>X), (11) 

where the divergence of X(f> was given by equation (R) above. Stokes Theorem and the relation m = dH A fi yields 
then 

H<f>= [ d{p<f>) = / du [ fj,div(X(f>). (12) 



lH=E JH<E J-oo JH 

Hence, taking a further energy derivative : 

d 



. )E IXe{4>) = M&v{<t>X)), (13) 



as we wanted to show. 



III. GENERAL REMARKS 



• Non-uniqueness: Note that the vector field X is far from unique. One may add any vector field tangent to the 
energy surfaces. This corresponds to a reparametrization of the surfaces and does not change the average of 
div (<pX). It does, however, affect fluctuations in this observable and hence a wise choice of X could lead to 
better convergence in numerical experiments. 

• Setting <f> = 1 it follows from (^J) and ((f); E) = 1 that : 

— — = (divX;E). (14) 



Neither this formula nor those in the Theorem make any reference to metric properties of M.. A metric 
on M allows one to compute gradients of functions and this gives one natural choice of the vector field X, 
namely X = VH/\\ VH || 2 . In the case of Euclidean space and Lebesgue measure this leads to the formula 
l/T{E) = (V • jp^p; E) as was found in @. 

• Iteration: Given the analytic expression for a vector field X and of <f> as in the Theorem the above formulas may 
be iterated indefinitely. Thus formula (|^) implies that any energy derivative of a ^-average can be measured 
within the /i-ensemble. In particular, any derivative of the temperature, e.g. specific heat, is also measurable in 
the /i-ensemble. 

• Fluctuations: Using the identity div((j>X) — {X ■ V)0 + (j) divX we may rewrite the equation (^) as follows : 

^-(0; E) - ((X ■ V)0; E) + (8<t> 5{dwX)- E). (15) 



where 5(f) = <p — {(f>; E) etc. Thus the energy derivative has a contribution coming from the derivative of (f> in the 
direction of X as well as from the product of fluctuations in the observables (f> and divA, cf. example B below. 
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IV. EXAMPLE A : THE CANONICAL ENSEMBLE 



As an illustration of the Theorem we will show that the canonical ensemble can be measured through u-averages. 
The canonical ensemble is most conveniently here defined as the Laplace transform of the /i-ensemble ||j4j • More 
precisely, the weighted partition function is given by : 



Z(j3,<t>)= \me- pH <i>= J dEe-^VHstt). (16) 

Assuming that ^_e(</>) has an analytic extension in E and s is small enough we may expand the partition function by 
iterating ((?]). It is convenient to write Dx(4>) = div (4>X) in terms of which 

LLE+M) = e s ^ f i E (cf > ) =VE(e sDx <f>)- (17) 
Here the last exponential should be interpreted through its power series. Therefore, formally : 

Z(f3,4>) = J dse-^ E+s ^ E (e sDx 0), (18) 

where in this last expression E has a fixed value. Dividing (normalizing) by the factor /xe(1) = e s ( E ' we see that all 
the terms in the expansion of the exponential can be measured in the /^-ensemble. Thus, up to a constant factor the 
canonical partition function can in principle be evaluated in the ^-ensemble. Of course, this approach does raise the 
question about convergence. In general one can only hope for the expansion to be asymptotic, e.g. in the number 
of degrees of freedom in the system. To make our statements more explicit consider the standard partition function 
(0=1) and the following expansion of the /^-entropy : 

S(E) = S(E ) + J2 {E ~ k f 0)k ^ ( 19 ) 

fc>0 

As a consequence of Theorem 1 all the numbers Sk (but not S(Eq)) are measurable in the /j,-canonical ensemble. We 
shall assume in the following that S2 is strictly negative^ The partition function we write as : 

Z(f3) = [dE e^ E+s ^ = e- F ^\ (20) 



Now, let f3 be close to si = dS / dE(E$) and let V(j3) be the extremal point (the conjugated variable in the Legendre 
transform) of the exponent —j3E + S(E). Again it is possible to calculate T(/3) perturbatively from the known 
quantities Sk and equation ([l9|). T(/3) satisfies (cf. also Bailyn || Section 11.6]) : 

=8E^ l -8E^W (21) 
Inserting E — T((3) + z we get for the exponent : 

- W) + S(T(P)) + 0(r(/3))z 2 i + V(z) (22) 

where V(z) = «3Z 3 /3! + «4Z 4 /4! + ... can also be calculated in terms of the coefficients Sfc. As the second derivative was 
assumed negative, we may carry out the Gaussian integral by standard techniques to obtain the (exact) asymptotic 
formula : 

F(J3) = F cl (P) - log[e^-w)^) 2 e^) lx ^ ], ( 23 ) 

with 



1 This corresponds to the assumption of thermodynamic stability, i.e. that the specific heat should be positive in the /j,- 
ensemble. At least in the thermodynamic limit there are mathematical arguments justifying such an assumption. 
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F cl (P) = const + (3T(f3) - S(T((3)) - ~ Iog(-^). (24) 

For the sake of clarity we shall in the following neglect the V^-term and use the last expression (corresponding to the 
classical action in QFT). Taking a p derivative we get for the average energy in the canonical ensemble : 

(H;0) = F^{0)=m-~^, (25) 
Differentiating once more we obtain the approximative formula for the specific heat : 

As a trivial but analytically accessible example (cf. jlO]]) we consider N harmonic oscillators for which the /^-entropy 
equals S{E) = (N - 1) logE and the extremum of -(3E + S(E) is attained for T(/3) = (N - l)//3. We get : 



N-l l-2(iV-l)//3 3 _ N 
1T~2 (7V-1)//? 2 ~J 



which happens to recover the exact canonical results. 



V. EXAMPLE B : INTERACTING PARTICLES 



Our second application is concerned with the thermodynamics of N particles in a box of volume V in R 3 . We shall 
compare the micro-canonical computations, e.g. of temperature, given here with formulae obtained from the canonical 
ensemble theory. It turns out that a natural assumption on the fluctuations in the kinetic energy is sufficient to obtain 
equivalence of the two approaches. 

The Hamiltonian is taken to be of the form 

3N 

H(p,q) = ^> 2 /2 + EM?) + UeM, (28) 
i=l 

with the standard symplectic structure on R 6N . Here K(p) = J2iPi/ 2 ^ s the kinetic energy, C/ cxt = 2j=i Uj(qj) 
is a box confining potential and U- m t is an interaction potential, e.g. a sum of two-body interactions. Note that for 
notational convenience the momentum vector is here considered as a 3./V dimensional vector whereas the configuration 
coordinates arc considered as N three-dimensional vectors. 



For the vector field X we choose : X\ = p/2K(p). The reader might worry about the fact that this vector field is 
singular at p — but as we shall see below the singularity is integrable when the number of particles is sufficiently 
large (at least 2 particles are needed in the applications below). A direct computation shows that : 

div (Xi) = (3N - 2)/(2K (p)). (29) 

The temperature is then given by : 

I 3N-2 



T(E) 



-{l/K{p)-E). (30) 



In order for the average to be well-defined it is necessary that 1/K(p) is integrable at p — 0, i.e. that J l/p 2 d 3N p < 00 
where the integral is over a neighborhood of p = 0. This is the case when 3iV > 2, i.e. when the system contains at 
least one particle. 

„:c„ u„„+ 1 /„ <9 TV en _ TV E^2 d 



The inverse specific heat, I/c = J^T(E) = — T(E) 2 Jfg yTe\ > can ^ e calculated using (^) from which 



y T{E) dE> T(E) v 1 K h 1 4 \ / > / v ; 
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and thus 

(3N-4)(l/K 2 :E) 

Vc(£) = l- {m _ 2)[l/K , E? - ( 32 ) 

This time we need that 1/K 2 (p) is integrable and this happens when 3N > 4, i.e. at least two particles are present. 
The standard definition of the (canonical) temperature is 

T c (E) = ^(K(p);E). (33) 

Had the average been in the canonical ensemble this would indeed have been the canonical temperature. We have the 
following formula for the ratio : 

3iV — 2 

T c (E)/T(E) = ^ r {K;E}{l/K;E). (34) 

We may also calculate the inverse specific heat this time using the canonical temperature. From (jfy and the obvious 
identity div(KXi) = 3N/2 we obtain : 

1 - = (35) 



c c (E) dE ^ ' T(E) 

In order to compare the above formulas we shall consider the fluctuations in the kinetic energy, defined by 

K = (K(p);E) +6K. (36) 

Our assumption in the following will be that (K) = (K(p); E) and (5K) 2 are both of order N. In particular, that for 
large enough TV the singularities are integrable and we have the expansion : 



11 SK (SK) 



2 



K (K) (K) 2 (K) 3 
Taking the average on both sides yields : 



+ o(N~ 2 ). (37) 



X K' (K) (K) 
and similarly by squaring before taking the average 

1 . 1 , ((SK) 2 ) 



<i> = 7^ + ^+^- 2 ) (38) 



(T^> = TTTTa + 3 ^ -T- + o(N~ s ) (39) 



K K 2 ' (K) 2 (K)i 
Using d38|)-(|39]) and retaining only terms to order iV _1 both ( |32"1) and (|3^) reduces to 



c(E) 3N (K; E) V ' 

We have here given a micro-canonical derivation of an expression which relates the specific heat to fluctuations in the 
kinetic energy. This was previously found by Lebowitz et al Q using an ingenious technique of inverting ensemble 
averages. Their method, however, relies on (unstated) analytic properties of the ensembles involved whereas the 
method presented here makes the assumptions explicit in terms of the fluctuations in the kinetic energy. We also 
note that our assumption on SK gives a sufficient condition for the equivalence not only of T(E) and T C (E) (the ratio 
( |34| ) differs from 1 by a term of order TV -1 ) but also for the derived expressions for the inverse specific heat (to order 
iV _1 ). Even close to a phase transition where l/c(E) tends to zero we would expect that our assumptions (|3^)-(|39|) 
are not violated. 

As a final application consider the ^(-canonical pressure exerted by the N particles on the walls of the container. We 
may define it as the average force per surface area of the container. The pressure, P, can then be calculated through 
the Virial Theorem (see e.g. Becker p. 98] or Abraham and Marsden 0, Example 3.7.32]). The j'th particle 
j = 1, . . . , N having coordinates qj is confined by the external potential Uj(qj). The force, dF — P dA, exerted on the 



G 



surface element dA is given by the (time-) average of the external forces in a small neighborhood SV of that surface 
element, i.e. P dl = (]T\ VUj : q 3 6 SV). 

This formula is vectorial but a scalar quantity is obtained by taking the scalar product with the coordinate of the 
volume element in question. Summing over the whole surface removes the restriction on the g^'th coordinate, thus 
yielding 

P fq-dA=(£q j -VU j ) (41) 
J j 

and finally, using Stokes Theorem we see that the left hand side equals 3PV. By the Virial Theorem the time average 
of (pd p — qd q )H vanishes (since it is a total time derivative). Hence 

PV/N=±(2K(p)-<t>;E), (42) 

with $ = J2j=i Qj • dU int (q)/dqj. 

We therefore obtain the following derivative : 

-^-(PV/N) = c(E)^-(PV/N) = 1 - ^1 A (43) 
dT y 1 ' v ' dE y 1 ' 3N dE^ ' 1 v ' 

Equation ( |l5| ) implies the exact formula : 

-^{PV/N —T) = c (E)^y^(S$ 5(1/K);E), (44) 

which under the assumption above on the kinetic energy fluctuations reduces to 

■^-(PV/N — T)= . 2c V %so (5$ 5K-.E) + oil). (45) 
dT \ I ! (2,NT(E)) 2 ^ ' ' v 1 v ' 

This again is a micro-canonical rederivation of a formula previously obtained by Lebowitz et al Q under the afore- 
mentioned analyticity assumption of ensemble inversion. 

I am grateful to Cristian Giardina, John Milnor, Ettore Aldrovandi and Dietmar Salamon among others who have 
helped me clarifying the ideas and calculations presented here. 
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